BIOENGINEERING AND BIOTECHNOLOGY 



ORIGINAL RESEARCH ARTICLE 

published: 03 December 2013 
doi: 10.3389/fbioe. 2013. 00018 



Effect of red blood cells on platelet activation and 
thrombus formation in tortuous arterioles 

Jennifer K. W. Chesnutt 12 and Hai-Chao Han 13 * 

' Cardiovascular Biomechanics Laboratory, Department of Mechanical Engineering, The University of Texas at San Antonio, San Antonio, TX, USA 

2 Department of Pathology, University of Texas Health Science Center at San Antonio, San Antonio, TX, USA 

3 Biomedical Engineering Program, UTSA-UTHSCSA, San Antonio, TX, USA 



Edited by: 

Tim Harrigan, Johns Hopkins 
University Applied Physics 
Laboratory, USA 

Reviewed by: 

Zhi-Yong Li, Southeast University, 
China 

Seungik Baek, Michigan State 
University, USA 

Haiyan Luo, Cooper Union, USA 
Yang Liu, The Hong Kong Polytechnic 
University, Hong Kong 

'Correspondence: 

Hai-Chao Han, Department of 
Mechanical Engineering, The 
University of Texas at San Antonio, 
One UTSA Circle, San Antonio, TX 
78249, USA 

e-mail: hchan@utsa.edu 



Thrombosis is a major contributor to cardiovascular disease, which can lead to myocar- 
dial infarction and stroke. Thrombosis may form in tortuous microvessels, which are often 
seen throughout the human body, but the microscale mechanisms and processes are not 
well understood. In straight vessels, the presence of red blood cells (RBCs) is known 
to push platelets toward walls, which may affect platelet aggregation and thrombus for- 
mation. However in tortuous vessels, the effects of RBC interactions with platelets in 
thrombosis are largely unknown. Accordingly, the objective of this work was to determine 
the physical effects of RBCs, platelet size, and vessel tortuosity on platelet activation and 
thrombus formation in tortuous arterioles. A discrete element computational model was 
used to simulate the transport, collision, adhesion, aggregation, and shear-induced platelet 
activation of hundreds of individual platelets and RBCs in thrombus formation in tortuous 
arterioles. Results showed that high shear stress near the inner sides of curved arteriole 
walls activated platelets to initiate thrombosis. RBCs initially promoted platelet activation, 
but then collisions of RBCs with mural thrombi reduced the amount of mural thrombus 
and the size of emboli. In the absence of RBCs, mural thrombus mass was smaller in a 
highly tortuous arteriole compared to a less tortuous arteriole. In the presence of RBCs 
however, mural thrombus mass was larger in the highly tortuous arteriole compared to the 
less tortuous arteriole. As well, smaller platelet size yielded less mural thrombus mass and 
smaller emboli, either with or without RBCs. This study shed light on microscopic interac- 
tions of RBCs and platelets in tortuous microvessels, which have implications in various 
pathologies associated with thrombosis and bleeding. 

Keywords: erythrocytes, shear-induced platelet activation, tortuosity, thrombosis, microvessels, mean platelet 
volume, computational simulation, discrete element model 



INTRODUCTION 

Tortuous microvessels are found throughout the human body, 
such as coronary vasculature ( Hutchins et al., 1978), and cerebral 
(Spangler et al, 1994; Brown et al., 2002), retinal (Sasongko et al, 
2012), and conjunctival (Owen et al, 2008) arteries. Microvascu- 
lar tortuosity alters blood flow to increase fluid shear stress that 
can activate platelets and induce thrombosis, even in the absence 
of vessel injury and hypercoagulability. For example, an in vivo 
study showed that thrombus formed due to high shear stress in 
rat venules that were made curved from an originally straight 
shape, but did not form in straight venules (Liu et al, 2008). As 
well, microvascular thrombi have been observed in humans in 
clinical, experimental, and autopsy settings (Gando, 2010). Con- 
sequently, the study of platelet activation and thrombus formation 
in tortuous microvessels is of clinical importance. 

In straight vessels in vivo, large numbers of red blood cells 
(RBCs) flow in the vessel pushing platelets to the wall, known 
as platelet margination as demonstrated in vitro (Aarts et al, 
1988). Thus, platelets are highly concentrated near the lumen walls. 
Responses of platelet margination and platelet adhesion rates to 
RBC collisions, hematocrit, and platelet size have been analyzed 



using mathematical models (Tokarev et al., 2011a,b). Also, the 
microscale processes of RBC and platelet interactions have been 
analyzed using computational simulations tracking individual 
cells. These simulations have been used effectively to examine the 
effects of RBC collisions, RBC aggregation, RBC deformability, 
RBC cytoplasm viscosity, platelet size and shape, platelet adhesion 
forces, and channel size on platelet margination, shear forces on 
platelets, and the process of thrombus formation (Miyazaki and 
Yamaguchi, 2003; AlMomani et al, 2008; Mori et al, 2008a; Ches- 
nutt and Marshall, 2009a; Kamada et al., 2012; Reasor et al, 2012). 
Though these previous models have elucidated effects of RBCs in 
straight vessels and channels, effects of RBCs in tortuous vessels 
have not been addressed. In non-cylindrical vessels, a few in vitro 
studies with blood flow through sudden expansions showed that 
hematocrit affected spatial concentrations of RBCs and platelet- 
sized particles (Zhao et al., 2008) and affected platelet adhesion 
(Karino and Goldsmith, 1984). Hence, RBCs also likely play a role 
in thrombosis in tortuous vessels. 

Most computational simulations that modeled thrombus for- 
mation by tracking individual platelets were performed in straight 
tubes or channels without RBCs (Miyazaki and Yamaguchi, 2003; 
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Pivkin et al, 2006; Filipovic et al, 2008b; Fogelson and Guy, 2008; 
Mori et al., 2008b; Kamada et al, 2010), and very few included 
RBCs (Mori et al, 2008a; Xu et al., 2009; Kamada et al., 2012). Few 
simulations have been accomplished in non-cylindrical geome- 
tries (e.g., stenosis or tubular expansion), and these studies were 
without RBCs (Filipovic et al, 2008a; Kamada et al, 201 1). In most 
of these models, platelets were activated due to an injured segment 
of the vessel wall. The effects of platelet activation by high shear 
stress were not studied. 

Therefore, we have previously simulated the microscale 
processes of thrombus formation in tortuous arterioles and 
venules based on shear-induced activation of individual platelets 
(Chesnutt and Han, 2011, 2013). However, these studies did not 
include RBCs in the simulations. 

Another physical factor that is present in conditions associated 
with many thrombotic and bleeding complications is platelet size, 
which is measured clinically as mean platelet volume (MPV). For 
normal healthy human subjects, an increase in MPV was shown to 
increase platelet aggregation in platelet rich plasma (Karpatkin, 
1978). Elevated MPV is observed in pathological conditions, 
including diabetes (Papanas et al, 2004), hypertrophic cardiomy- 
opathy (Cambronero et al., 2009), acute myocardial infarction 
(Chu et al., 2010), restenosis following coronary angioplasty (Chu 
et al., 2010), pulmonary hypertension (Guvenc et al, 2012), and 
giant platelet disorders (Mhawech and Saleem, 2000). MPV also 
decreases under other pathological conditions, such as reactive 
systemic amyloid A amyloidosis (Erdem et al.,2012) and Wiskott- 
Aldrich syndrome (Ochs et al, 1980). Hence, it is also of interest 
to determine the effects of platelet size in the presence of RBCs. 

Consequently, the objective of this work was to determine the 
physical effects of RBCs, MPV, and vessel tortuosity on shear- 
induced platelet activation and thrombus formation in tortuous 
arterioles. These results would help to develop new treatment 
strategies and diagnostic tests for the risk of cardiovascular disease. 

MATERIALS AND METHODS 

The transport, collision, activation, adhesion, and aggregation 
processes of hundreds of individual RBCs and platelets of differ- 
ent MPV were numerically simulated in tortuous arterioles. The 
computational simulation conditions are described first, followed 
by brief descriptions of the employed discrete element model and 
shear-induced platelet activation model. 

COMPUTATIONAL SIMULATION CONDITIONS 

Segments of tortuous arterioles, each having diameter D = 25 \im, 
were modeled as two-dimensional (2D) channels in the shapes of 
cosine curves with two periods. Walls near the bends (locations 
with highest curvature) were designated as either inner or outer 
walls (Figure 1). A tortuosity index T was defined as the ratio of 
amplitude A to wavelength \ of the cosine curve (i.e., T = A/\). To 
determine the effect of tortuosity, a low tortuosity arteriole with 
T = 0.09 and a high tortuosity arteriole with T = 0.21 were exam- 
ined. Each arteriole had the same arc length. The assumption of 
2D flow was used for computational efficiency and was expected 
to give qualitatively similar results and similar relative effects due 
to RBCs, tortuosity, and MPV compared with three-dimensional 
(3D) flow. The chosen arteriole diameter, tortuosity indices, and 




FIGURE 1 | Schematic of the high tortuosity arteriole showing the 
locations of inner walls, outer walls, and bends, which are labeled 
numerically. 



shapes were similar to those observed in tortuous arterioles in 
humans (e.g., Brown et al., 2002; Sasongko et al., 2012). 

Platelets and RBCs were modeled as 3D spheres for compu- 
tational efficiency, similar to previous computational studies of 
thrombus formation that utilized spheres or circles (Mori et al, 
2008a; Kamada et al., 2010; Flamm et al, 2011). This approxi- 
mation was expected to have a small effect on the results and is 
justified in more detail in the Section "Discussion." In a given 
simulation, platelets had uniform properties, and RBCs had dif- 
ferent uniform properties. To determine the effect of platelet size 
in pathological states, simulations used three different values of 
MPV found in humans (3.8, 7.1, 15.8 fl). The volume of an RBC 
was set to the physiological value in humans of 94 fl. The densities 
of platelets and RBCs are similar, and so were set to the same cell 
density p c = 1063 kgm~ 3 . 

To determine the effects of both platelet distribution and pres- 
ence of RBCs, three different seeding scenarios were employed 
to seed cells (i.e., platelets and RBCs) into the arteriole. In all 
seeding scenarios, cells entered the arteriole at the inlet with an 
initial velocity equal to the local fluid velocity at the location of 
the cell centroid, and seeded platelets were initially unactivated. 
In the whole-lumen-seeded platelets scenario, initial positions of 
platelet centroids at the inlet followed a pseudorandom proba- 
bility distribution that excluded the region roughly one platelet 
radius from the wall such that platelets entered anywhere along 
the inlet boundary without touching the wall. RBCs were not 
included in the whole-lumen-seeded platelets scenario. In the 
near-wall-seeded platelets scenario, initial positions of platelet 
centroids at the inlet were bounded to regions within 7.5 |xm 
of the walls. Centroid positions within these regions followed a 
pseudorandom probability distribution, without platelets touch- 
ing the wall at their initial seed locations. This seeding assumed 
that the platelet distribution at the inlet was that expected in 
a straight arteriole with whole blood. RBCs, though, were not 
included in the near-wall-seeded platelets scenario. The RBCs 
with platelets scenario was the same as the near-wall-seeded 
platelets scenario, except that RBCs were also seeded into the 
flow. Positions of RBC centroids were initialized about the cen- 
ter of the arteriole, which spanned the 10-|xm-region in which 
platelets were not seeded. Centroid positions of RBCs within 
this region followed a pseudorandom probability distribution. 
In each seeding scenario after a cell was seeded at the entrance, 
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the cell was subjected to no other restrictions on its location or 
motion. 

The number of unactivated platelets that entered the arteriole 
per unit time (seeding rate) was set such that a physiological time- 
averaged platelet count (300,000 mm~ 3 ) would be achieved in the 
absence of platelet adhesion to the walls. Hence, the seeding rate 
for whole-lumen-seeded platelets was required to be about twice 
that for near- wall-seeded platelets and for RBCs with platelets. The 
physiological time-averaged hematocrit was set to 6%. This hema- 
tocrit is representative of that in microvessels, which can be much 
lower than the ~40% hematocrit of arteries and veins (Keller et al, 
1994; Lipowsky,2005). 

Simulations were performed for Newtonian, incompressible, 
steady, fully developed, 2D flow in the x-y plane. Steady flow 
was assumed because flow pulsatility is minimal in the micro- 
circulation. Due to low values of hematocrit and volume fraction 
of platelets, the effects of cells, aggregates, and mural thrombi 
on the fluid flow were neglected. This assumption is reasonable 
since the focus of this work was the initial formation of throm- 
bus, rather than later possible occlusion. Centroids of 3D spherical 
cells remained the x-y plane in the 2D flow because forces on cells 
in the z-direction were zero, and the only non-zero torques on 
cells were in the z-direction. A Poiseuille velocity profile with typ- 
ical mean arteriole velocity [7 = 0.6cms _1 was imposed on the 
fluid at the inlet of the arteriole. In a straight arteriole, this veloc- 
ity profile would yield a wall shear rate lyo ' = 1440 s -1 , which 
is in the range of typical values in arterioles (Kroll et al., 1996). 
Simulations considered blood to be a suspension of RBCs and 
platelets in plasma. Considered as a single-phase fluid, whole 
blood is a non-Newtonian fluid with variable apparent viscos- 
ity that depends on shear rate, which is mainly due the presence 
and aggregation of RBCs (Skalak et al, 1981). Rather than treat 
blood as a continuous fluid, the computational model simulated 
both plasma and cells, and incorporated a model that accounted 
for increased drag force on a cell due to increased apparent viscos- 
ity from cell crowding in the vicinity of the cell (Di Felice, 1994). 
Because of these modeling assumptions and the fact that plasma 
is a Newtonian fluid (Chien et al, 1966), the fluid was chosen 
to be Newtonian plasma with fluid density pf = 1030 kg m -3 and 
dynamic viscosity u, = 1.2 cR The parameter values utilized gave a 
Reynolds number (Re = pf UD/\i) ~0.1 and particle Stokes num- 
bers (St= p c d 2 !7/18|xD) on the order of 10~ 5 for platelets and 
10~ 4 for RBCs, where d is diameter of either a platelet or RBC. 
Simulation parameters related to physical properties of arterioles, 
cells, and plasma, as described above, are listed in Table 1. 

Measures of thrombosis based on platelet activation, adhesion, 
and aggregation were determined throughout the simulation time. 
An aggregate was defined as a group of two or more platelets in 
contact with each other. A platelet was considered to be part of a 
mural thrombus if the platelet was in contact with the wall or if 
it was part of an aggregate that was in contact with the wall. One 
measure of thrombosis was mural thrombus mass, which was the 
number of platelets in mural thrombi at a given time multiplied by 
the mass of a platelet. This measure represented platelets that were 
in the vessel at the specified time. Another measure was average 
embolus mass, which was the number of platelets in an aggregate 
that was not in contact with the wall, averaged over the simulation 



Table 1 | Simulation parameters related to physical properties of 
arterioles, cells, and plasma. 



Parameter 


Value 


Comments 


D 


25 |xm 




1 (lOvv) 


n no 


Low tortuosity arteriole, 
/\^16|im, \^177|xm 


T (high) 


0.21 


High tortuosity arteriole, 
a ~ ou |xm, a. ~ 141 (jim 


MPV (small platelet) 


3.8fl 


Observed in Wiskott-Aldrich 
syndrome, corresponding to 

l.ci^lIM S|JMcl ICdl UldlllclfcM 


MPV (normal platelet) 


7.1 fl 


Physiological value, 
corresponding to 2.4 p.m 
spherical diameter 


MPV (large platelet) 


1 5.8 f I 


Observed in diabetes, 
corresponding to 3.1 |im 
spherical diameter 


Volume of an RBC 


94 f I 


Corresponding to 5.6 |im 
spherical diameter 


Platelet count 


300,000 mm -3 


0.2% Platelet volume fraction 


Platelet seeding rate 
(who e-lumen-seeded) 


0.67 ms" 1 


- 


Plate et seeding rate 
(near-wall-seeded) 


u. J4 ms 




Hematocrit 


D /o 


\.l 1 nib nDL btifcJUIliy 1 ale 


Pc 


I UDo Ky 111 




Pf 


lUoi) Kg m 






\.Z Cr 




u 
IyoI 


0.6 cms -1 
1440 s- 1 




Ycrit 


1500s~ 1 




Re 


0.1 




St (platelet) 


O (10 -5 ) 




St (RBC) 


O(10" 4 ) 





time and multiplied by the mass of a platelet. Activation count 
was defined as the cumulative number of activation events. This 
measure was a count of the number of platelets that had become 
activated since the start of the simulation. To compare seeding sce- 
narios that had different platelet seeding rates, activation percent 
was defined as the percent of seeded platelets that had become 
activated since the start of the simulation. 

Simulations were run to a final time ( tf = 1 .67 s) at which mural 
thrombus mass reached a state of statistical equilibrium for cases 
with the RBCs with platelets scenario. Simulation cases are listed in 
Table 2. Each simulation was performed as a serial (non-parallel) 
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Table 2 | Simulation cases. 



Case 


Tortuosity 


Platelet size 


Seeding scenario 


1 


Low 


Norma 


Whole-lumen-seeded platelets 


2 


Low 


Norma 


Near-wall-seeded platelets 


3 


Low 


Normal 


RBCs with platelets 


4 


High 


Normal 


Whole-lumen-seeded platelets 


5 


High 


Normal 


Near-wall-seeded platelets 


6 


High 


Norma 


RBCs with platelets 


7 


High 


Small 


Whole-lumen-seeded platelets 


8 


High 


Small 


Near-wall-seeded platelets 


9 


High 


Small 


RBCs with platelets 


10 


High 


Large 


Whole-lumen-seeded platelets 


11 


High 


Large 


Near-wall-seeded platelets 


12 


High 


Large 


RBCs with platelets 



program on a 3.0-GHz node having either four or eight processing 
cores, with 2 GB of RAM per core. Simulations were run simul- 
taneously on these processing cores. CPU times ranged from 8 to 
17 days per simulation. 

DISCRETE ELEMENT COMPUTATIONAL MODEL 

The transport, collision, activation, adhesion, and aggregation of 
hundreds of individual platelets and RBCs were simulated in tortu- 
ous arterioles by a mesoscale, discrete element method for blood 
cells. Details of the computational model have been previously 
published (Chesnutt and Marshall, 2009a,b; Chesnutt and Han, 
201 1). A brief description of the model is provided as follows. 

The computational model followed the motion of each spher- 
ical cell due to interactions with the fluid, vessel walls, and other 
cells. Activated platelets were subjected to adhesion with the wall 
and aggregation with each other. RBCs were subjected to aggre- 
gation with each other, though shear rates were high enough in 
this work to suggest disaggregation of RBC aggregates (Schmid- 
Schonbein et al, 1968). In the remainder of this subsection, the 
term adhesion is used to refer to both adhesion and aggregation. 

Cells were modeled as 3D spheres, rather than 2D circles, so 
that they possessed mass and could be subjected to the dominant 
3D forces and torques as defined below. Fluid flow was modeled in 
2D to increase computational efficiency. Cell centroids were ini- 
tialized on the 2D flow x-y plane. As the fluid velocity component 
in the z-direction was zero, and the only non-zero fluid vorticity 
component was in the z-direction, centroids of 3D cells remained 
in the x-y plane. Thus, 3D angular and linear momentum equa- 
tions were reduced to 2D and were solved for each 3D cell in 2D 
flow, given by 



dv dQ. z 
m— = F F + F A , I z —7- 
at at 



M F>Z + M A)Z , 



(1) 



where mis cell mass, vis cell velocity, 7 Z = (l/W)md 2 is moment of 
inertia, Q. z is cell rotation rate, Fp and Mp ;Z are respectively fluid- 
induced force and torque on the cell, Fa and Ma ;Z are respectively 
force and torque on the cell due to adhesion and collision, and t 
is time. Fluid-induced force Fp was composed of the dominant 
fluid-induced forces, which were drag force and added mass force. 



Adhesion and collision forces acting on two colliding cells con- 
sisted of a force due to elastic deformation of cells F ne , a force due 
to energy dissipation from collision F n( j, and a resistance force 
F s due to sliding of one cell over another. Adhesion and colli- 
sion torques acting on two colliding cells consisted of a torque 
generated by the sliding resistance force F s and a torque due to 
resistance of rolling of one cell over another M r . Adhesion and 
collision forces and torques were given by 



F A = F ne n + F nd n + F s t s , M A)Z 



d 



r F s (n x t s ) + M r (f r x n), 

(2) 



where n is a unit vector tangent to the line connecting the centroid 
of the cell to the centroid of the other cell, t s is a unit vector in 
the direction of relative motion of the cell surfaces projected onto 
the plane orthogonal to «, and t r is a unit vector in the direction 
of rolling velocity of the cell. Forces (drag force, added mass force, 
F n d, F s ) and torques (M FjZ , M r ) were determined using equations 
given in the literature. Elastic deformation force F ne was deter- 
mined, as described below, based on the theory of Johnson et al. 
(1971). 

Upon collision, cells were assumed to retain their shapes except 
within a flattened circular region of contact, in accord with the 
model of cell-cell contact of Bell et al. (1984). Receptor-ligand 
bonds that formed within the contact region were modeled as 
springs following (Bell et al., 1984). Under our model assumptions, 
receptor-ligand binding took a form that was mathematically anal- 
ogous to van der Waals adhesion with a time-dependent adhesive 
surface energy density a(t) given by 



a(t) = --(Xb 



Xe) / N b (t- t 0 (s)) sds, 



(3) 



where a is spring constant, x b is length of the bond, x e is equi- 
librium length of the bond, N b is number density of bonds, and 
s is radial position with the contact region divided by the radius 
of the contact region. The variation of N b with time was given by 
a kinetics equation of Bell (1978) that accounted for forward and 
reverse reactions of bond formation given by 



dN b 
dt 



k { (N L0 - N h ) (Nro - N h ) - krN b , 



(4) 



where kf and k r are respectively forward and reverse reaction rate 
coefficients, and JVlo and Nro are respectively initial ligand and 
receptor densities. The reaction rate coefficients were given by 

Dembo et al. (1988) 



k { = kg) exp 
k x = /c r o exp 



2kf 



(a - Ots) (*b - Xe) 

2kt 



(5) 



where kf 0 and fc r o are respectively initial forward and reverse equi- 
librium reaction rates, a ts is transition state spring constant, k is 
the Boltzmann constant, and T is absolute temperature. Number 
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density of bonds in Eq. 4 was obtained using the solution given 
by Chesnutt and Marshall (2009a). For simulations in the current 
work, the time scale for bond formation was much less than the 
time scale for elastic response of two colliding cells. Thus, num- 
ber density of bonds was approximated by its equilibrium value 
Nb(oo), such that adhesive surface energy density simplified to a 
constant value, given by 



a(f) 



(Xb - ^) 2 N b (oo), 



(6) 



where calculations yielded a(t) = — 4 x 10 _1 dyncm -1 for 
platelets and — 1 x 10~ 4 dyn cm -1 for RBCs. 

Because of the analogy with van der Waals adhesion, we used 
the theory of Johnson et al. (1971) for adhesive elastic particles 
to approximate elastic deformation forces acting between two 
colliding cells, given by 



^ne 

On 

Sc 



4 («o) 4 («o) 

~(-) 2 -*-(-) 



(7) 



where a(t) is contact region radius, and 8 n is normal overlap dis- 
tance of cell surfaces. Critical force F c and critical overlap S c can 
be written in terms of contact region radius at equilibrium «o as 



2(6 1 /3)fl' 



\ E 



(8) 



where _R and E are respectively effective radius and effective elastic 
modulus of the two cells, with adhesive surface energy density a 
given by Eq. 6. 

Mechanisms of platelet aggregation and adhesion differ 
depending on shear rate (Jackson et al, 2009). As maximum 
shear rates in the current work ranged from about 1 800-3 1 00 s~ 1 , 
platelets were assumed to adhere to von Willebrand factor immo- 
bilized on the subendothelium through glycoproteins GPIb or 
GPIIb-IIIa, and to aggregate with each other through bridging of 
GPIb or GPIIb-IIIa by plasma von Willebrand factor (Ikeda et al, 
1991; Kroll et al, 1996; Konstantopoulos et al, 1997; Shankaran 
et al, 2003). At low shear rates, RBCs aggregate with each other, 
while at high shear rates these aggregates break up. Though 
different theories of the mechanisms of RBC aggregation exist 
(e.g., Chien and Jan, 1973; Neu and Meiselman, 2002), models 
of receptor-ligand binding of RBCs used in computational stud- 
ies have produced aggregation behavior similar to experimental 
observations (Bagchi et al., 2005; Chesnutt and Marshall, 2009a, 
2010). Hence, in the current work, RBC aggregation was assumed 
to occur through receptor-ligand binding. 

Typical ranges of values obtained from the literature for para- 
meters in Eqs 4 and 5, which vary widely for different cell adhesion 
problems, are given in Table 3, along with values used in simula- 
tions. Few values of parameters in Table 3 are known specifically 
for platelet adhesion or aggregation. We selected the density of 
platelet receptors and other parameters to achieve an effective 



Table 3 | Values of parameters related to receptor-ligand binding used 
in simulations and ranges of experimental values. 



Parameter (units) 


Value used in 


Value used in 


Range of 




simulations for 


simulations 


experimental 




platelets 


for RBCs 


values 


i i 9 1 i 

fcfo (cm' s 1 ) 


10 ' 


10 ' 


10" 12 to 10" 7a 


k, 0 (s" 1 ) 


10" 5 


10" 5 


10" 5 to 10 b 


W L0 , W R0 (cm" 2 ) 


1.115 x 10 12 


4.798 x 10 9 


10 9 to 10 l2c 


x e (nm) 


20 


20 


5-50 d 


x b (nm) 


29 


29 


0 to length at 








which bond breaks 


a (dyn cm -1 ) 


1.98 


2.2 


0.01-10 e 


cr ts (dyn cm -1 ) 


-5 


-5 


-5 to 5 f 


t(K) 


310 


310 





"From Hammer and Apte 11992). 

"From Bell (1978), Ward and Hammer (1993). 

'From Demboetal. (1988), Doggettetal. (2002), Fteiningeretal. (2006), Michelson 
(2007), Fogelson and Guy (2008). 

"From Bell et al. (1984), Evans and Leung (1984), Dembo et al. (1988). 
"From Dembo et al. (1988), Chtcheglova et al. (2004). 
'From Dembo et al. (1988). 

spring constant for adhesion and aggregation of a platelet that 
was of the same order of magnitude as that determined in a 
previous study. Specifically, the effective spring constant for a 
normal sized platelet at the equilibrium contact region radius, 
calculated as rjTta 2 Nb(oo), was 34Nm~'. This value compared 
well with the platelet effective spring constant 50 Nm -1 , which 
was determined by Filipovic et al. (2008b) from computational 
simulations by comparison with experiments. Parameters for 
RBCs were chosen to obtain the constant adhesive surface energy 
density a = — 1 x 10~ 4 dyn cm -1 , which previously showed good 
agreement with in vitro aggregation experiments (Chesnutt and 
Marshall, 2009a). 

SHEAR-INDUCED PLATELET ACTIVATION MODEL 

We used our model of shear-induced platelet activation as in 
Chesnutt and Han (201 1), which was shown in simulations to ini- 
tiate thrombosis at inner walls in agreement with previous in vivo 
experiments. This activation model ensured the assumption that 
a platelet cannot become activated by a physiological shear rate or 
shear stress, irrespective of the amount of time a platelet is sub- 
jected to physiological shear stress. As well, this activation model 
clearly indicated the point at which an individual platelet became 
activated. These two attributes of our activation model are not 
necessarily apparent in other existing models that also account for 
exposure time to shear stress. 

Our model assumed that a platelet became activated if it expe- 
rienced a shear rate above a critical shear rate (Ycrit) defined as 
Ycrit=/lYol, where / is a scaling factor greater than unity. We 
assumed a critical shear stress of t = 60 dyn cm~ 2 to determine 
critical shear rate. With approximate apparent viscosity of whole 
blood (ib = 4 cP, the critical shear stress corresponded to a critical 
shear rate y cr i t = x/|Xb = 1500 s _1 . Hence, the scaling factor was 
chosen to be/= 1.042. The critical shear stress was chosen based 
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FIGURE 2 | Shear rate contours of the high tortuosity arteriole (top) 
and low tortuosity arteriole (bottom). 



A B 




6" TOO 200 300 (Mm) 0 TOO 201) 1pm) 



FIGURE 3 | Thrombus formation at t, with different seeding scenarios, 
including whole-lumen-seeded platelets (top), near-wall-seeded 
platelets (middle), and RBCs with platelets (bottom), for low tortuosity 
index (A) and high tortuosity index (B), for Cases 1-6 Cells include 
activated platelets (black), unactivated platelets (white), and RBCs (red). 
Platelets are normal size. 



on in vitro experiments with whole blood and platelet rich plasma, 
which showed that shear-induced platelet activation and aggrega- 
tion occurred at shear stresses that ranged from 15 to 30 dyn cm~ 2 
(Chow et al., 1992; Konstantopoulos et al., 1995, 1997; Shankaran 
et al., 2003) . Although some of these experiments also showed large 
increases in platelet activation and aggregation as shear stresses 
surpassed 75-140 dyn cm -2 , our critical shear stress was chosen 
to illustrate trends in the early stage of thrombus formation with 
respect to changes in shear stress due to tortuosity. 

To account for the presence of chemical agonists released by 
activated platelets, the model assumed a platelet became activated 
if it contacted another activated platelet, similar to the model of 
Kamada et al. (2010). In our simulations, only activated platelets 
were subjected to adhesion and aggregation, and activation was 
considered irreversible. 

RESULTS 

EFFECTS OF RBCs AND TORTUOSITY WITH PLATELETS OF NORMAL SIZE 

In both arterioles, maximum shear rates exceeded the wall shear 
rate that would occur in a straight arteriole (lyol = 1440 s _1 ) and 
the critical shear rate for activation (lycritl = 1500 s _1 ). Maximum 
shear rates were 1764 and 3099 s _1 in the low and high tortuosity 
arterioles, respectively. The highest shear rates occurred at inner 
walls in both arterioles (Figure 2), and fluid velocities in these 
regions were higher for the high tortuosity arteriole (not shown). 

Cases with platelets of normal size for the two tortuosity indices 
and three seeding scenarios (Cases 1-6) are presented in this sub- 
section. In each case, platelets first became activated at locations 
where shear rates were greater than the critical shear rate (critical 
shear regions) of the first and second bends. The first few mural 
thrombi formed on the inner walls of these two bends in some 
cases (Cases 2, 4, 6), and in other cases formed downstream of 
these two bends before reaching the third bend (Cases 1, 3, 5). 
After the onset of a thrombus at inner walls of the first or second 
bend, the thrombus first grew along the wall in the downstream 
direction and then later in the upstream direction. Throughout the 
simulations in each case, mural thrombi continued to form along 
the walls at other locations, as shown at the final time in Figure 3. 
For cases with RBCs, almost immediately after the first few mural 
thrombi formed, RBC collisions with the thrombi caused thrombi 
to roll along the wall. Throughout the simulations with RBCs, 
RBC collisions caused thrombi to continually roll along the wall 
and exit, generally in a single layer of platelets that could con- 
tain 10 or more platelets. For cases without RBCs, thrombi did 
not roll along the walls, except for one activated platelet and one 
three-platelet mural thrombus in Case 5. 

For each case, emboli continually detached as mural thrombi 
grew into the flow (away from the wall), with some emboli reat- 
taching on walls downstream. A single layer of platelets remained 
attached to the wall just after an embolus detached from a mural 
thrombus. For cases without RBCs, emboli consisted of an aver- 
age of 8.4-8.9 platelets, with a maximum of 49 platelets. However, 
with RBCs, emboli were smaller and consisted of an average of 2.6 
platelets, with a maximum of 6 platelets. 

Activation count initially became larger for cases with RBCs 
compared with cases with only near-wall-seeded platelets, for a 
given tortuosity index (Figures 4A,B). This promotion of platelet 
activation due to RBCs lasted for a longer time in the high 



tortuosity arteriole (~0.6 s) than in the low tortuosity arteriole 
(~0.2 s). However, after this short time, activation count became 
smaller with RBCs compared to near-wall-seeded platelets with- 
out RBCs. Due to the higher seeding rate with whole-lumen- 
seeded platelets, activation count for this seeding scenario became 
larger than the near- wall-seeded platelets scenario and the RBCs 
with platelets scenario for both tortuous arterioles, as shown in 
Figure 4C for the high tortuosity arteriole. The low tortuosity 
arteriole showed similar results to those in Figure 4C, except 
that activation counts at the final time were slightly smaller at 
623, 517, and 462, with whole-lumen-seeded platelets, near- wall- 
seeded platelets, and RBCs with platelets, respectively. These same 
values for the high tortuosity arteriole were 654, 532, and 490. 
However, activation percent was lowest at most times with whole- 
lumen-seeded platelets (Figure 5), due to the seeding of some 
platelets near the center of the arteriole. 

For a given tortuosity index at any given time after 0.23 s at 
most, mural thrombus mass was smallest with RBCs compared 
to cases without RBCs, and about the same compared between 
cases without RBCs (Figure 6). Toward the end of the simulations, 
mural thrombus mass was generally larger for the low (versus high) 
tortuosity arteriole in the presence of RBCs, but generally smaller 
for the low tortuosity arteriole in the absence of RBCs. 
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FIGURE 4 |Time variation of activation count during an early stage, for Seeding scenarios include whole-lumen-seeded platelets (dashed-dotted 
the low tortuosity (A) and high tortuosity (B) arterioles; and during the line), near-wall-seeded platelets (dashed line), and RBCs with platelets (solid 
entire simulation time for the high tortuosity arteriole (C), for Cases 2-6 line). Platelets are normal size. 
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FIGURE 5 | Activation percent versus tortuosity index, with 
whole-lumen-seeded platelets (triangles), near-wall-seeded platelets 
(squares), and RBCs with platelets (circles), at an early time t = 0.2 s 
(A) and the final time t, = 1.67 s (B), for Cases 1-6. Platelets are normal 
size. 
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FIGURE 6 | Time variation of mural thrombus mass, with 
whole-lumen-seeded platelets (blue line), near-wall-seeded platelets 
(black line), and RBCs with platelets (red line), for the low tortuosity 
arteriole (A) and high tortuosity arteriole (B), for Cases 1-6. Platelets are 
normal size. 



EFFECTS OF RBCs AT DIFFERENT PLATELET SIZES 

Cases with platelets of different sizes (MPVs) for the three seed- 
ing scenarios in the high tortuosity arteriole (Cases 4-12) are 
presented in this subsection. In these cases, the first few mural 
thrombi formed either at the inner walls of the first and second 
bends or on walls downstream of bends after platelet activation 
at the first and second bends. The location of formation of the 
first few mural thrombi did not have a clear relationship with 
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FIGURE 7 | Thrombus formation at f = 1.54sfor different seeding 
scenarios, including whole-lumen-seeded platelets [top), 
near-wall-seeded platelets {middle), and RBCs with platelets [bottom), 
for small platelets (A) and large platelets (B), for the high tortuosity 
arteriole for Cases 7-12. Cells include activated platelets (black), 
unactivated platelets {white), and red blood cells (red). 



MPV or seeding scenario. Similar to the cases of normal MPV, 
RBC collisions with abnormal MPV caused mural thrombi to roll 
along walls and exit the arteriole. Additionally, without RBCs, most 
cases with abnormal MPV exhibited a small amount of rolling of 
thrombi, with the large MPV case exhibiting some rolling of large 
groups of platelets, though less than cases with RBCs. Locations of 
cells at the end of the simulations are shown for abnormal MPV 
in Figure 7. 

Emboli continually detached from mural thrombi for each case, 
with some emboli reattaching on walls downstream. Emboli in the 
absence of RBCs consisted of an average of 6.6-9.5 platelets (max- 
imum of 37 platelets), and in the presence of RBCs an average of 
2.3-2.8 platelets (maximum of 7 platelets) (Figure 8A) . For a given 
seeding scenario, although the average number of platelets per 
embolus decreased with increasing MPV, average embolus mass 
increased with increasing MPV (Figure 8). 

With near-wall-seeded platelets of large MPV, mural thrombi at 
the first bend eventually grew upstream along the wall and reached 
the inlet boundary at t = 1.54 s. In this case, locations of new near- 
wall-seeded platelets could not be seeded into the flow, and so 
simulations were analyzed at the time when thrombi first blocked 
the inlet near the wall, at t = 1.54 s, instead of the final simulation 
time t = 1.67 s. In future studies, the computational model can be 
modified to seed near-wall-seeded platelets toward the center of 
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the arteriole in the event that locations near the walls are occu- 
pied by thrombi. Nevertheless, at time t = 1.54 s, mural thrombus 
mass had reached a state of statistical equilibrium for cases with 
the RBCs with platelets scenario so that results were expected to 
be valid. In cases without RBCs, mural thrombi at inner walls of 
the first and second bends grew in both the upstream and down- 
stream directions. With RBCs, the small MPV case showed slight 
upstream growth of mural thrombi at the first and second bends 
before thrombi began to roll along the wall, while the large MPV 
case showed no upstream growth before thrombi began to roll 
along the wall. 

Unlike the case with normal MPV near the beginning of the 
simulation, activation count with abnormal MPV was not clearly 
larger with RBCs compared to only near-wall-seeded platelets. 
As simulations progressed, activation count was largest for the 
case with whole-lumen-seeded platelets, followed by the case with 
near-wall-seeded platelets, and smallest for the case with RBCs, 
for a given MPV (e.g., Figure 9A). Activation percent at the end 
of the simulation was smallest for the case with whole-lumen- 
seeded platelets due to possible seeding in the center of the 
arteriole (Figure 9B). At the final time, activation count (and acti- 
vation percent) increased as MPV increased for a given seeding 
scenario. 



For a given MPV, mural thrombus mass increased with time 
with similar values for the two cases without RBCs and was small- 
est with RBCs. The trends with abnormal MPV (not shown) were 
similar to those with normal MPV (Figure 6B), except that values 
with small MPV were generally smaller and values with large MPV 
were generally larger, at a given time for a given seeding scenario. 
For example, values at the final time were at most 53% lower with 
small MPV (0.3-0.8 ng) and at most 176% higher with large MPV 
(1.0-3.6 ng). Though at the final time for a given seeding scenario, 
for cases with abnormal MPV compared with normal MPV, the 
number of platelets in mural thrombi did not vary as much as 
mural thrombus mass nor it have as clear a relationship with MPV 
(Figure 9C). Final numbers of platelets in mural thrombi with 
abnormal platelets varied from 23% lower to 25% higher than 
with normal platelets, depending on MPV and seeding scenario. 

DISCUSSION 

RBCs, MPV, AND TORTUOSITY AFFECT THROMBUS FORMATION 

Current results suggested that for normal sized platelets at the 
onset of thrombosis, RBCs collided with unactivated platelets to 
push them closer to walls where critical shear regions or mural 
thrombi were located. After a very short time, this effect was likely 
mitigated by formation of large enough mural thrombi that unac- 
tivated platelets could easily contact. Later, RBC collisions caused 
mural thrombi to roll out of the arteriole for all platelet sizes, 
which reduced platelet activation events, mural thrombus mass, 
and embolus sizes. Collisions of RBCs with mural thrombi caused 
emboli to detach before mural thrombi could grow as large as in the 
case without RBCs. However, without RBCs, mural thrombi were 
able to grow larger toward the vessel center where larger emboli 
were detached from mural thrombi by larger fluid forces, in com- 
bination with collisions with other platelets. The overall effect of 
RBCs on thrombus formation in tortuous arterioles may be to 
reduce the occurrence of thrombotic occlusion of the arteriole or 
of vessels downstream. 

Results of the two seeding scenarios without RBCs indicated 
that thrombosis (and hemostasis in the case of injury) would 
proceed in a similar manner, but would be more efficient with 
near-wall-seeded platelets in a tortuous arteriole, as also expected 
in a straight vessel. Because the near-wall-seeded platelets sce- 
nario virtually excluded collisions of RBCs within the arteriole 
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FIGURE 8 | Average number of platelets per embolus versus MPV 
(A) and average embolus mass versus MPV (B), with whole-lumen- 
seeded platelets {triangles), near-wall-seeded platelets {squares), and 
RBCs with platelets {circles), for the high tortuosity arteriole for Cases 
4-12. Values were averaged over the entire simulation time. 



800 



600 



o 400 



=8 200 




5 10 15 
MPV (fl) 



10 15 
MPV (fl) 



10 15 
MPV (fl) 



FIGURE 9 | Activation count versus MPV (A), activation percent versus MPV (B), and number of platelets in mural thrombi versus MPV (C), with 
whole-lumen-seeded platelets {triangles), near-wall-seeded platelets {squares), and RBCs with platelets {circles), in the high tortuosity arteriole at 
t= 1.54 s for Cases 4-12 



Frontiers in Bioengineering and Biotechnology | Biomechanics 



December 2013 | Volume 1 | Article 18 | 8 



Chesnutt and Han 



RBCs affect thrombosis in arterioles 



while it mimicked the effect of RBCs on the concentration pro- 
file of platelets at the inlet, a partial effect of RBCs would be to 
increase the efficiency of thrombosis and hemostasis in tortuous 
arterioles. However, because the RBCs with platelets scenario (i.e., 
near-wall-seeded platelets scenario with the effect of RBC colli- 
sions) caused less mural thrombus mass and smaller emboli sizes, 
our results suggested that RBCs would increase the efficiency of 
thrombus initiation, but reduce the risk of thrombotic occlusion 
due to arteriole tortuosity. 

Our work suggested that, due to their small diameter, platelets 
of small MPV initially had a lower probability of collision with 
RBCs that could push them toward walls, and so RBCs did not 
increase activation count initially, as in the case of normal MPV. 
Moreover, RBCs did not increase activation count initially in the 
case of large MPV, which was possibly due to platelets of large 
mass undergoing smaller rebound velocities after collisions with 
RBCs that could push them toward walls. Our study showed that 
increased MPV yielded increased activation count at the final time 
for a given seeding scenario, including the two scenarios without 
RBCs and the scenario with RBCs. This result agreed with our 
previous study with whole-lumen-seeded platelets without RBCs 
(Chesnutt and Han, 2013), which was demonstrated at the tortu- 
osity indices used in the current study, in addition to two inter- 
mediate tortuosity indices. Larger MPV may promote thrombosis 
through larger mural thrombus mass and larger emboli that would 
be more likely to occlude the arteriole or vessels downstream. 

Our results revealed that tortuosity could generate high enough 
shear rates to induce platelet activation. Our results also suggested 
that, in the absence of RBCs, mural thrombi remained more intact 
in the arteriole of low tortuosity than in the arteriole of high tortu- 
osity, likely due to forces at inner walls of bends (fluid and collision 
forces) that were lower in the arteriole of low tortuosity. Therefore, 
the fact that mural thrombus mass was smaller in the presence of 
RBCs in the arteriole of low tortuosity must be attributed to one or 
more of the following presumed events: more RBC collisions with 
mural thrombi, higher RBC collision forces with mural thrombi, 
or more rolling of mural thrombi out of the vessel at low tor- 
tuosity. Additional simulations to specifically measure individual 
collision events and collision forces or to virtually exclude specific 
events, such as thrombus rolling, may elucidate reasons for these 
results in future studies. 

Our simulation of the process of thrombus formation yielded 
results that were similar to previous in vivo, in vitro, and computa- 
tional studies. For example, thrombus formation was initiated at 
inner walls in some cases. This result agreed with a previous in vivo 
study in rats, which showed thrombosis was initiated at inner walls 
of curved venules due to high shear stress and high shear stress gra- 
dient in those regions (Liu et al, 2008) . Another study with porcine 
whole blood perfused through glass stenoses that were coated with 
collagen showed that platelets initially adhered at the downstream 
end of the stenosis with some adhering upstream of the apex of the 
stenosis (Para et al., 201 1). Later, thrombus formed both upstream 
and downstream of the stenosis apex. Although the very high shear 
rates in this previous study (~ 1 00,000 s -1 ) indicated that platelet 
adhesion and aggregation occurred through different mechanisms 
than in our study (see Materials and Methods), the process of ini- 
tial thrombus formation was similar to our study, in that thrombi 



that formed at the region of highest shear grew downstream first 
and then upstream later. 

In a previous 2D computational study in simple shear flow with 
shear rate of 1000 s~\ platelets became activated by an injured wall 
segment to form a thrombus (Mori et al., 2008a). RBCs caused the 
thrombus to grow more along the wall than into the flow, and 
caused larger a thrombus mass. Our simulations agreed with the 
former result but not the latter. The discrepancy may be partially 
due to differences in the width of the region in which platelets 
were seeded, which was larger in the previous study (17 \x,m) than 
in our study (7.5 u,m). Our smaller seeding width provided less 
space for thrombus to grow away from the wall and more oppor- 
tunities for RBCs to collide with the layer of platelets adhered to 
the wall, causing thrombus to roll and exit the vessel. However, 
in the previous study of Mori et al. (2008a), RBCs were farther 
away from the thrombus base such that the thrombus was several 
platelet layers deep, and RBCs collided with platelets to push them 
toward the thrombus, rather than detach them from the wall. 

A recent 3D computational study simulated individual platelets 
and deformable RBCs in straight cylindrical tubes of diameters 
between about 12 and 15 |xm, with length 40 u,m and Re about 
0.01-0.02 (Kamada et al, 2012), which were smaller vessels and 
flow rates than in our study. In these previous simulations with 
and without RBCs, platelets became activated due to an injured 
wall segment, adhered to the injured wall, and initiated throm- 
bus formation. The mural thrombus then grew downstream, with 
some platelets detaching from the thrombus due to collisions with 
RBCs or due to fluid forces. These initial processes of thrombus 
formation and effects of RBC collisions and fluid forces were in 
agreement with our results, though thrombosis was initiated by 
high shear stress in our study, rather than injury. However, the 
presence or absence of RBCs did not change thrombus height, area 
of injured wall covered by the thrombus, or number of platelets 
in the thrombus in this previous study. As noted by the authors, 
this result was due to the deformation of RBCs during collision 
with the thrombus. In contrast, our study assumed rigid RBCs 
and found smaller mural thrombus mass in the presence of RBCs. 
However, our Reynolds number (Re = 0.1) was 10 times higher 
than in the previous study, such that if deformable RBCs were 
used in our study, we would have observed larger RBC velocities 
that might produce large enough RBC collision forces to detach 
platelets from mural thrombi. Differences in results between the 
presence and absence of RBCs would be less pronounced with 
deformable RBCs than with rigid RBCs, but we would expect an 
effect of RBCs. 

MODEL LIMITATIONS 

The model of shear-induced platelet activation did not account 
for exposure time to shear stress, which has been shown to affect 
activation in vitro (Heliums, 1994). Our previous work (Chesnutt 
and Han, 2013) showed that our activation model was either in 
quantitative or qualitative agreement with a shear- induced platelet 
activation model that accounted for magnitude of and exposure 
time to shear stress (Alemu and Bluestein, 2007), depending on 
the parameters chosen for the model of Alemu and Bluestein. 

The model contained various assumptions for computational 
efficiency. Platelets and RBCs were assumed to be spherical and 
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rigid in the absence of collision or adhesion. Although unacti- 
vated platelets are discoid, activated platelets are more spheroidal 
(Michelson, 2007). Hence, the shapes of activated platelets are 
closer to spheres than disks. As the focus of our work was on 
thrombus formation by activated platelets, the difference in results 
due to use of a spherical platelet shape, rather than a spheroidal 
shape, were expected to be small. RBCs have a resting shape that is a 
biconcave disk, and can exhibit either a tumbling motion or tank- 
treading motion in fluid flow depending on shear rate (Bitbol, 
1986). As well, RBCs are highly deformable, which is necessary for 
RBCs to traverse capillaries, which have smaller diameters than 
RBCs (Skalak et al., 1989). Because our study was simulated in 
arterioles with diameters larger than capillaries and focused on 
the effect of the overall presence of RBCs on thrombus formation, 
the assumption of spherical non-deformable RBCs was expected 
to have a small effect on the results. 

Fluid flow was assumed to be 2D and unaffected by the presence 
of cells (one-way coupled). Under these assumptions, our previ- 
ous work showed that aggregates of cells in 2D shear flows were 
quantitatively similar in size and shape to those in 3D flows (Ches- 
nutt and Marshall, 2010), and that thrombosis was initiated at the 
same locations in 2D channel flows as in vivo (Chesnutt and Han, 
2011). Also, the effects of cells on the fluid were reduced due to 
the low volume fractions of platelets and RBCs in our simulations. 
Because this work focused on initial thrombus formation, rather 
than occlusion, the assumptions of 2D flow and one-way coupling 
were expected to give qualitatively similar results and similar rel- 
ative effects due to RBCs, MPV, and tortuosity compared to 3D 
flow and two-way coupling. 

CLINICAL RELEVANCE 

The physical effects of RBCs, MPV, and tortuosity may compound 
the chemical factors and other physical factors present in dis- 
ease states. In the absence of RBCs, our results showed much 
larger mural thrombus mass and embolus size, which suggested a 
higher risk of thrombotic complications. This result was in accord 
with a retrospective study of anemic patients who underwent free 
flap reconstructive surgery (Hill et al., 2012). Preoperative val- 
ues of hematocrit and hemoglobin were significantly lower in 
patients who experienced thrombosis and who experienced flap 
failure, which was primarily due to thrombosis. The authors pro- 
posed that thrombosis could have been partly due to increased 
turbulence from lower blood viscosity or due to alterations in 
platelet aggregation and concentrations of coagulation factors. In 
addition to these factors, our results suggest that fewer RBC colli- 
sions due to lower hematocrit might promote thrombosis once 
thrombi begin to form during free flap procedures in anemic 
patients. 

In cases with RBCs, mural thrombus mass and embolus mass 
were about double for cases with large platelets compared to cases 
with normal platelets. Hence, the physical effects of large platelets 
and microvascular tortuosity in diabetes may be partly responsible 
for the large incidence of thrombosis, which is also attributed to 
hypercoagulability in diabetes. Because metformin reduces platelet 
size in diabetes (Dolasik et al, 2013), this drug may also reduce the 
risk of thrombosis due to the physical effects of platelet size. This 
mechanism was speculated in our previous work that neglected 



RBCs (Chesnutt and Han, 2013) and is confirmed by the current 
work that includes RBCs. Mural thrombus mass was about half for 
cases with small platelets compared to cases with normal platelets, 
throughout the simulations with RBCs. Assuming a similar effect 
under conditions of platelet activation due to vessel injury, rather 
than high shear stress, this physical effect of small platelet size 
may contribute to the high risk of bleeding observed in Wiskott- 
Aldrich syndrome, which is also attributed to severely low platelet 
count and other factors (Ochs et al., 1980). 

CONCLUSION AND OUTLOOK 

The computational model simulated individual platelets and RBCs 
in the microscale processes of shear-mediated thrombus forma- 
tion in tortuous arterioles. We concluded that RBCs initially pro- 
moted platelet activation, but that this effect was quickly reversed 
as simulations progressed. Collisions of RBCs with mural thrombi 
yielded reduced mural thrombus mass and embolus mass. As well, 
a reduction in platelet size resulted in reduced mural thrombus 
mass and embolus mass, either with or without RBCs. Hence, ves- 
sel occlusion may be less likely due to the presence of RBCs and 
due to smaller platelet size. A smaller mural thrombus mass was 
observed at the low tortuosity index compared to the high tor- 
tuosity index in the presence of RBCs. Our results suggest that 
the physical mechanisms associated with RBCs, platelet size, and 
vessel tortuosity play important roles in platelet activation and 
thrombus formation. 

To our knowledge, except for our two previous works (Chesnutt 
and Han, 201 1, 2013), no other computational studies have simu- 
lated shear- induced activation of individual platelets in thrombus 
formation, possibly because previous activation models did not 
clearly indicate the incident of activation. With the clear occur- 
rence of platelet activation in our computational model, we were 
able to simulate a large number of platelets in thrombus formation. 
Our study shed light on mechanical factors involved in tortuosity- 
induced thrombosis, including the presence of RBCs, platelet size, 
and microvessel tortuosity. The model serves as a foundation to 
build upon in future work, for example with improvements such 
as more physiological shapes of RBCs and platelets and/or 3D fluid 
flows. Future studies may help determine the risk of thrombotic 
or bleeding complications due to physical mechanisms associated 
with RBCs, platelets, and vessel tortuosity, to aid treatment and 
prevention of these complications in various pathologies. 
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